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QQ ' It was recently proposed to use small groups of trapped ions as qubit carriers in miniaturized electrode ar- 

^SJ , rays that comprise a large number of individual trapping zones, between which ions could be moved 1 1 2^. 

This approach might be scalable for quantum information processing with a large numbers of qubits. Pro- 
cessing of quantum information is achieved by transporting ions to and from separate memory and qubit 
^ ' manipulation zones in between quantum logic operations. The transport of ion groups in this scheme plays 

, a major role and requires precise experimental control and fast transport. In this paper we introduce a theo- 

. retical framework to study ion transport in external potentials that might be created by typical miniaturized 

' Paul trap electrode arrays. In particular we discuss the relationship between classical and quantum descrip- 

' tions of the transport and study the energy transfer to the oscillatory motion during near-adiabatic transport. 

, Based on our findings we suggest a numerical method to find electrode potentials as a function of time 

, to optimize the local potential an ion experiences during transport. We demonstrate this method for one 

' specific electrode geometry that should closely represent the situation encountered in realistic trap arrays. 
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^ . 1 Introduction 

Quantum information processing is a rapidly evolving field of physical science. Its practical importance 
arises from the exponential speedup in computation of certain algorithmic tasks over classical computation 
|3 1. Building an actual device that can process quantum information, however, is technologically difficult 
■ due to the need for qubits that can be processed and read out with high fidelities and the extreme sensitivity 

of the quantum mechanical states stored in these units against external uncontrolled perturbations. A 
promising technical approach as shown over the last decade, is to use strings of ions as physical qubits 
confined in Unear electromagnetic Paul traps LL^. These strings are stored in a single trap and constitute 
a one dimensional crystallized structure whose vibrational modes can be laser cooled to their ground states. 
The strong mutual coupling of the ions by Coulomb forces in such a crystal has been proposed and utilized 
to create arbitrary superpositions of quantum states of the ionic internal states {[^5_,6\). In the last 
few years methods were developed that enable quantum state engineering with high precision and long 
coherence times fT T, TTO* TTj . The necessary criteria fTT\ for large-scale quantum computation have 
been demonstrated in the past years, and small algorithms have been implemented successfully [13. 14, 15l 
1161 1171 . However, as in other approaches aiming towards quantum computation, scaling to many qubits 
is challenging. Considerable overhead is required by quantum error correcting schemes that permit robust 
quantum computation and make large-scale implementations feasible. To scale up a linear string of many 
ions, a rapidly growing number of vibrational degrees of freedom needs to be controlled and cooled to the 
ground state for reliable processing. This is extremely difficult to realize. A more recent proposal p||2) 
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has been made to circumvent this problem by using small arrays of a few qubits that are shuttled around in 
two-dimensional microstructures to process and store quantum states at various locations. 

An initial systematic study showed that coherent transport of ions in linear trap arrays is possible with 
nearly no loss in contrast during the motion 1 18 1. In this experiment an adiabatic transport of a qubit was 
performed over a distance of 1.2 mm in a time span of about 54 /is with negligible heating. CuiTently, there 
are strong efforts under way to demonstrate the possibility of building large-scale ion trap structures. For 
example, suggestions have been made to combine miniaturized ion chips directly with CMOS electronics 
to handle the resources required to control the many electric potentials L19J . Moreover, fast transport 
requires excellent experimental control of all these potentials. 

A detailed scheme of how a viable architechture of an ion trap processor could look has been recently 
studied by Steane |20|, fully incorporating quantum error coiTecting codes. The physical gate rate of this 
proposed 300 qubit processor unit was found to be limited by 

2 10 

Tg = 1 h Tcool + Tp a few /iS, (1) 

with the first two terms being an average time of the part of a typical gate that involve motion that is times 
for splitting 2i'^qj^), recombining and moving lOi^jT^) a small ion string, where t'coM and i^^ are 
typical axial and radial trapping frequencies, respectively. The last two terms correspond to cooling after 
the transport has been done, and the time duration of conducting the actual phase gate, respectively. On the 
other hand, if large amounts of energy are transferred to the ions, longer cooling times might be needed. 
Inserting typical operating conditions shows that the first two terms make up a considerable part of the 
performance of the physical gate rate. In order to keep this part as small as possible we need designs for 
electrode structures enabling fast qubit transport. 

In the following we present a theoretical framework that governs the transport dynamics of ions trapped 
in a time varying external potential. In Sect. 2 the equations of motion for the transport are derived. Sect. 
3 discusses the general classical solution in terms of an Ermakov parametrization. This approach is useful 
to express the quantum approach presented in Sect. 5, which uses the Heisenberg picture following the 
approach of Kim et al. 1211 . In Sect. 4 we point out some well-known properties of a quantum harmonic 
oscillator exposed to a transporting force for the simpler case when its frequency is kept constant. Sect. 
5 presents the general quantum solutions and the interrelation between classical and quantum transport. 
Based on this framework we discuss in Sect. 6 a well-controlled regime for the transport and also include 
first order perturbations to the transport dynamics. In Sect. 7, we present numerical optimization routines to 
extract optimum switching of potentials for the transport and study miniaturization of electrode structures 
to estimate the required resources for a well-controlled transport. Finally, a simple electrode model is used 
to find a practical rule for the segmentation of ion traps revealing insight into the resources needed for 
large-scale layouts, that should be also apphcable for more general trap arrays. 



2 Classical equations of motion 

A linear segmented Paul trap, e.g. as used in recent experiments P lll4l [131116112211231 . consists typically of 
two alumina wafers with gold coated electrode surfaces of a few micrometer thickness. The slotted wafers 
provide electrical RF and DC fields for 3D confinement of ions. The arrangement for control electrodes is 
schematically sketched in Figlsjwhere only a single layer is shown. The confinement along the a;-axis is 
achieved solely by electrostatic fields whereas the remaining two orthogonal radial directions correspond 
to a dynamical trapping by ponderomotive RF forces. In this article we limit ourselves to transport along a 
single dimension x from —6/2 to 6/2. If we denote the coordinate of the ion in the laboratory frame by q 
then we have from Newton's equation of motion 

m + qi-t,)=' qi-to)^0, (2) 

m oq 2 
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with two initial conditions as the equations on the rhs; Q is the elementary charge and m the mass of the 
transported ion. We assume a time interval and location of the ion starting at —to and —6/2, and ending 
at +to and 6/2, respectively. In order to make use of coherent states of a harmonic oscillator (that do not 
spread in time) we are interested in designing the time-dependent electrical potential as 

(j){q,t) - 0iocai(9 - qait)) + </?o(i), (3) 

where (t>\o<:-i\{q) ~ muj^q^ /2Q is purely quadratic with constant curvature in a sufficiently large range 
around the minimum, and f{){t) is a time-dependent offset with no influence on the dynamics. Here, 
we prescribe the dynamics by specifying a desired harmonic frequency ojq and the temporal shift of the 
harmonic well by a transport function qQ{t). The residual, uncontrolled force caused by insufficient flex- 
ibility in creating the desired harmonic potential deteriorates the transport performance. Its effect can be 
described by the difference potential or residual acceleration, i.e. 

(/-les <?, t) = (j){q, t) - ipo{t) and a,^,{q) = , (4) 

2Q m oq 

respectively. Due to imperfect realization of the harmonic well (pre^iq, t) adds fluctuating parts to the ideal 
harmonic potential as a function of position or time, critically depending on the electrode structure used. 
In section 7 we will discuss a numerical scheme for approximating (l){q, t) based on superpositions of 
individual electrode potentials in an optimal way. 

We finally can write down the classical equation of motion 

il + ujqU = ~qo{t) + a,s^{u + qait)) (5) 

which we transformed into a frame moving with qQ{t) by u = q ~ qg. The net acceleration on the rhs 
corresponds to an external force and displaces the ion from its equilibrium position u = in this frame. 
Since we will treat only the first two perturbation terms we expand the final equation of motion around the 
minimum of the well and rearrange some terms to get 

il + u;^(l- a'^Jqo{t)]/ujl^u - a'JJqo{t)]/2 + . . . = -qo{t) + a,,,[qo{t)] (6) 

where primes denote differentiation with respect to u. For the following discussion we abbreviate uj^{t) = 
uJq{1 — a(.gJ(7o(i)]/^o) write f{t) = —qo{t) + ares[qo{t)] for the rhs of Eq.(|6|l. For certain electrode 
structures, we can disregard terms involving the second and higher order derivatives of a^esiq) (cf. Section 
7). We will make this assumption throughout the paper. In that case Eq.(|6j simplifies to the equation of 
motion of a parametrically driven and forced harmonic oscillator with the Hamiltonian 

^(t)^f + (7) 
2m 2 



and p = rait. 



3 Classical dynamics of ion transport 

To obtain a general classical solution with an arbitrary frequency modulation we first consider the for- 
malism which is most often used in conjunction with time-dependent invariants within so called Lewis- 
Riesenfeld methods |30|. These approaches have been shown to be successful in the quantization of time- 
dependent harmonic oscillators with many different kinds of time-dependencies. Here, we discuss the gen- 
eral classical solution using the Ermakov equation and its generalized phase equation for time-dependent 
frequencies. We then employ in section 5 the approach of Kim et al. 121J to express the general quantum 
solution in terms of its classical solution. 
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3.1 Homogeneous Solution 

Neglecting higher order terms we find the homogeneous part of the solution of Eq.(|6|l by setting /(t) = 0, 
thus solving 

Uc + w2(t)Mc = (8) 

for an arbitrary time-dependent frequency Lu{t). For this, it is most convenient to make the ansatz 

= p(i)e*A'(*) U2 ^ p{t)e-'^^'\ (9) 

introducing an amplitude function p{t) and a phase function /^(t), both real. Inserting Eq.(|9j into Eq.(|8} 
and considering real and imaginary parts results in the two equations 

p- pfi^ +uj'^{t)p^Q, 2pii + pfl = Q. (10) 

p is an integrating factor for the second equation on the right so that we can write 

= (11) 

where we have chosen the integration constant as 1. The constant on the rhs of Eq. JI 1> has the SI units 
rad/s that should be taken into account at the end. If we substitute this back into the first equation of 
Eq.lQO} we obtain the Ermakov equation for the amplitude function p{t) 

p + u;\t)p^l/p\ (12) 

For periods of constant frequency w = the general solution is ' 

p{t) = ±uj^^''^ v^cosh 6 + sinh 5 sin(2woi + d) , (13) 

where 5, 9 are constants of integration, their values depend on the past evolution 1301 . The solution for the 
generalized phase is easily obtained once p is known. From Eq. lll U we have 

p{t) = /* dt'p{t')-\ (14) 

J -to 

The general homogeneous solution is then given by 

Uh{t) ^ acp{t) cos{p,{t) + if), (15) 
with the classical amplitude Oc and initial phase if fixed by the initial conditions. 

3.2 Green's function and general solution 

We use the general framework of Green's functions to define a particular solution to the inhomogeneous 
case of Eq.(|6}, where we again terminate the expansion, i.e. a[2^ [q\ = for n > 2, to stay in a harmonic 
regime. Using the two independent homogeneous solutions of Eq.(|9j we can determine the causal Green's 
function 

G{t, t') = e(t - t')p{t)p{t') sm{p{t) - fi{t')), (16) 
with 9{t — t') the Heaviside function. Employing G{t, t'), a particular solution is given by 

up{t)^ f dt'G{t,t')fit')^p{t) f dt'sini^l{t)-p{t'))pit')f{t'). (17) 

J-to J-to 



' The general solution of this equation is easily obtained by first using p as an integrating factor with the integration constant 
2uj cosh <5. This equation is immediately transformed to a harmonic oscillator by x = — coshS/uj. 
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For later convenience we define the auxiliary function = ie J c?i'e*^'* ^ p{t')f{t') which will 
be useful for expressing the general quantum solution. In this notation we can abbreviate the particular 
solution by Up{t) = p{t) {C{t) + C*(t)} /2. 

Thus, we obtained the general solution as the sum of the general homogeneous solution Eg. dlSt and a 
particular solution 

u,(t)=Uh(t)+Up(i) = pW/2{ace*^''^*^+'^)+CW}+c.c. . (18) 

Higher derivatives, like velocity and acceleration, can easily be found from the general solution in Ea.(ll8> 
by using the Leibniz rule. For initial conditions where we start in the classical ground state Oc = we 
define the quantity S(i) = Up{t) + iuj{t)up{t), assuming that the transport starts at —to > —oo, i.e. later 
than the infinite past, and demand that it takes a finite amount of time. With the help of the last definition 
the energy transferred to the oscillator at instants ti (where /(ti) = 0) is then given by 

W(ti)=m|S(ii)|V2, (19) 



with 

/ (itV(<')(p(ti)/i(<i)cos(A^itO + {p(ii)+«^(ti)p(ti)}sin(A/iitO)/(i'), (20) 

J —to 

and Apit' = p{ti) — p{t'). We will call 5(t) the adiabatic suppression amplitude and its absolute square 
the adiabatic suppression factor Thus, we have derived the classical energy transfer for arbitrary frequency 
evolutions and arbitrary external transport forces. To evaluate this expression one first must solve for the 
explicit time-dependence of p and p, according to Eqs.( I12I1'41 i by integrating the Ermakov equation and the 
phase equation, and finally compute the transferred energy at different times using EQS.( ll9l2'0t . 



3.3 Adiabatic limit 

Since we are mainly interested in an adiabatic solution we can simplify the last expression by considering 
adiabatic expansions of the homogeneous solution for a parametrically driven harmonic oscillator. We 
introduce an adiabatic time scale T such that 



dj/uj - for T > 1 



(21) 



The general adiabatic expansion of the differential equations Eqs. J12ll4> is readily obtained 1291 
1 1 3 



Pit) 



and p{t)^uj{t)--^:^ + -^r^T + 



(22) 



This procedure is equivalent to a perturbative approach on the first term in Eq. (ll2t 1301 . We require that 
at instants —to,ti, i.e. at times when we measure the oscillator's energy, the frequency has settled into 
a constant. Also for the following discussion we define that the oscillator's initial frequency at —to is 
uj{—to) = ujQ, so that pi—to) = Xj ^JZjo, p{—to) = 0. Taking into account only the lowest order of the 
expansion in Ea.(l22> the expression in Ea.(l20> reduces to 

E{ti) = ^/uj{ti) f ' dt'f{t')uj{t')-^^^e'^''^^' with Apw = / ' dr uj{t) (23) 

J-to Jt' 

providing the adiabatic energy transfer in the first order of frequency modulation. 
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4 Quantum and classical, dragged harmonic oscillators with constant fre- 
quency 

Husimi [271 and Kerner f281 independently considered the forced quantum mechanical oscillator and found 
exact analytical expressions for their wavefunctions and propagators. We review some of their early ideas 
because they provide insight into the close relationship of the quantum and classical solution. In this 
paragraph we assume the frequency is independent of time. The corresponding Hamiltonian is given by 
Eq.0 with uj{t) — ujq. 

Following Husimi and Kerner, we can "uncouple" the classical oscillation by the transformation 

t) = 4>{u', t) exp{im,Ucu' /h), (24) 

with u' — u — Uc and Uc at first undefined. Inserting Ea.( l24> into the time-dependent Schrodinger equation 
gives 



^ (~^^^ ^™^o"'^^ (l)+m{uc+ujQUc-f)u'(f>-{rn/2){ul~Ljlul + 2fuc)(f>. (25) 



On the rhs we see that we can make the second term vanishing if we choose Uc to satisfy 

Uc + UJqUc - / = 0, 

i.e. if Uc satisfies the classical solution of Eq.Q. With this choice one can easily identify the classical 
action L{t) = {m/2){u^ — oj'^u^ + 2fuc) of a forced harmonic oscillator in the third term on the rhs of 
Ea.( l25> . Furthermore, if we make the ansatz 



dt'L{t') 



we can absorb this term as a time-dependent phase into 0. The remaining part of the wavefunction, %, 
then needs only to obey the usual harmonic oscillator wave equation in the frame defined by the classical 
trajectory with its internal coordinate u' 



In this way one can achieve a separation of the forced harmonic oscillator from the unforced oscillator in 
a frame moving with the classial trajectory. The wavepacket does not become deformed by the homoge- 
neously acting force. The quantum solution becomes displaced and only a phase is accumulated. 

To determine further properties we can assume now a stationary state with energy e„ for the solution of 
Eq.(|26j 

Xn{u', t) = Un{u') exp(-ie„t) e„ = (n + \)nwa, 

and evaluate transition probabilities at time t for the oscillator to be in the number state if it was initially 
in the number state m„ 



2 



Husimi and Kerner showed that these transition moments can be evaluated analytically 

Pmn{t) = (Ai!/i^!)7""^e-^(i;;-''(7))2 with 7(0 = m/2fi^o W + iuJUc? (27) 
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by using generating functions for the Hermite polynomials |'27"281. In Ea.(l27>. v is the greater while /i 
is the lesser of m and n, respectively. L^^~^ denote the associated Laguerre polynomials, and its time- 
dependent argument 7(t) describes the classical energy transfer in units of Thoq. From Ea. (l27> we see the 
classical character of the quantum solution: the transition probabilities are solely defined by the classical 
quantity ^{t). Also, if we consider starting from the ground state n = jj. = Q and using io(7) = 1 the 
probability distribution P„,o becomes a poissonian, and thus we find the signature of a coherent state. 

With this relation the expectation values for the mean energy and the dispersion of the energy distribu- 
tion are then immediately obtained 

{E^)n = fiwo ( ^ TO P„„ + 1/2) = fiu;o + 7 + 1/2) £« + ^01 (28) 

{{AE.mf)n = {hujof ((m - (to))2) = {nuJof{2n + 1)7 = 2e„ fi^o7 , 

where e„ is the initial energy before the force acts on the wavepacket. This is indicated in Ea.( l28> by the 
subscripts on the Ihs. Corresponding expressions for the classical solution 

{E)e„ =Eo + W {iAEf)Eo = '^EoW (29) 

are found if we average over the initial classical phase that are completely analogous to the quantum 
solutions.^ Eq is the classical energy before the transport and W = hujQ^{t) the classical energy transfer. 
The mean energy and the energy spread increase linearly with the energy transfer in both solutions although 
the energy distributions of the classical and quantum solution are quite different |27 1. Also, the zero point 
energy makes a difference between the classical and quantum description. If the system is initialized in its 
quantum ground state, transport can create a dispersion of the wavepacket due to eq > 0, while this is not 
the case for the classical ground state, i.e. if £^0 — 0. 



5 The dragged quantum harmonic oscillator 

Many methods have been developed to find exact quantum states of time- dependent oscillators. The gen- 
eralized invariant method by Lewis and Riesenfeld 1301 has been very successful in finding exact quantum 
motion in terms of wavefunctions and propagators. For the interpretation of time-dependent quantum sys- 
tems and for showing its relationships to their classical solution, however, the Heisenberg picture is more 
appropriate since the Heisenberg operators for position and momentum obey similar equations of motion 
than the corresponding classical quantities. In this paragraph we aim to interpret the quantum solution 
using its classical analogue and therefore use the general approach of Kim et al. |21 1 that is based on the 
general invariant theory but acts in a Heisenberg picture, in contrast to the original approach. 

The general invariant theory starts out by defining an invariant operator I{t) that satisfies the Heisenberg 
equation of motion. Ji et al. |24 1 used a Lie algebra approach to find the most general form of the solution 
with some integration constants c^, i — 1,2,3, arbitrary defining the initial conditions (see discussion at 
the end of this paragraph). If we fix these parameters according to the conditions of Eq.(3.4) in f251 the 
generalized invariant is of the form 

lT{t)=uJi(B\t)B{t) + ^, (30) 

with LOj as a constant of motion, and the annihilation and creation operators are 

^ \/y { ^p^' - ^^^^ ~ ^} + '^^^^^ 
B\t) = ^[ (p-i + ^p) m ~c}- ^-^m- (31) 



^ This result is easily derived by averaging the general classical solution in Eg. II 81 given in Sect. 4 for a constant frequency over 
the phase interval [0, 2n]. 
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It holds that [B{t), (t)] = 1, where B{t), B^{t) are solely represented by the classical quantities p = 
p{t),iJ. — fJ,{t),( = ({t) as introduced in previous paragraphs. q{t),p{t) refer here to the Heisenberg 
operators for position and momentum and we have assumed in addition that /(—to) 0. 

From the Heisenberg equations of motion for B, i.e. dB{t)/dt — —i[B{t),T-C{t)], one can obtain the 
simple time evolution for these annihilation and creation operators 

B{t) = e-'^''■*^B{^ta) B^t) = e'^''^*^ B\~to) , (32) 

with the phase function. Their evolution in time is a simple time-dependent phase-shift mediated 
by the generalized classical phase referenced to the initial time —t^. This last property guarantees the 
time-independence of the invariant and the equivalence to the Hamiltonian (if /(— ^o) — 0) at the time 
-to: 

lT{t) = lTi-to)=ni-to). (33) 

Following Kim et al. we can equate hermitian and anti-hermitian parts on both sides of Eas. (l32> by using 
the relations Eas. (l31> to determine the time-dependent Heisenberg operators for position and momentum 



qit) = p{t) I q{-to)V^cos^iit) + iLJ^ sin fi{t) } + Up{t) (34) 



p{t) ^q{-to)ni^/uj^ [p{t) cos fi{t) - p{t) ^ sin p.{t)\ 

+ [p(t)-i cosp{t) + p{t) sin^(t)l + mup{t), (35) 



where (7(— to),p(— to) denote position and momentum operator at time —t^, respectively. Similar as Up{t), 
the classical velocity can be expressed as 2up{t) = {{p — ip~^)C + {p + *P^^)C* }■ 

Our chosen initial conditions, Ea.(l33t. cast momentum and position operators into their standard form 



'<-'"> = ^f^ + ^'> Pi-V'—^^iB-B'] (36, 

taking h = 1. Kim et al. define a more general Fock state space based on number states of the invariant 
rather than on the Fock state space of the Hamiltonian and point out its importance and advantegeous 
properties. However, due to our choice of the initial conditions these two state spaces are identical and 
their distinction is irrelevant for our discussion. We can define the Fock basis in the usual way by taking 
the operators at —to according to 

B^^ ftl 

|n)s = |n,-to)B with \n,t)B = ^=^\0,t)B, (37) 



where the vacuum state \0,t)B is extracted from _B(t) |0)b = 0. Furthermore, we introduce the time- 
independent coherent states in this Fock basis 



— \/ n ' 



n— ^ 



with the complex amplitude a = jale"*"^, because these states are the closest quantum equivalent to the 
classical solution and include the oscillator ground state for a = {). With these definitions the expectation 
values for the Heisenberg position and momentum operators from Eg s. ( 1341331 can be calculated using 
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Eq.(|36) and Eq.lEl 



{a\qit)\a) = J-p{t)\a\ cos[fi{t) + ^] + Up{t) 

V TO 

{a\p{t)\a) = \/2^|a||p(t)cos[^(t) + (^] ~ p-\t) sm[^l{t) + ip]j + mitpit) (39) 
= rn^{a\q{t)\a). 

This way we retrieve exactly the same form for the mean values of position and mometum for the quantum 
solution as we obtained in Eq.(ll8> for the classical solution. If we disregard the zero point energy in 
{a\T-C{—to)\a) / LUf) — jap + 1/2 ~ |ap and set the matrix element equal to the potential energy at a 
classical turning point, we have Oc ~ V2a;o|a| making the homogeneous solution of the classical and 
quantum formulations and hence the total solution identical. Here, xq = y/h/uom is the extension of 
the ground state wave function of the harmonic oscillator. Alternatively, a full quantum description in the 
Schrodinger picture can be obtained by employing the time evolution operator that can be represented as a 
product of time-dependent displacement and squeezing operators 1 2 1 27 1 . 
Similarly we can compute the dispersions of q{t),p{t) in the coherent state 

{a\{Aq{ma) = /2m {a\{Ap(t)f\a) = {p-' + p')m/2. (40) 

From the dispersion for the momentum we see that the wavepacket generally spreads solely due to the 
presence of the terms p^ and fp. These matrix elements do not depend on the force because the force acts 
homogeneously in space and equally on the whole wavepacket. After periods of frequency modulations 
the dispersions in Ea. (l40> are both time-dependent and exhibit oscillatory behaviour revealing a certain 
amount of squeezing |24|. For example, if we assume that after the transport we end up with a nonzero 
S we can use the exact solution in Eq. (I13> and evaluate the rhss of Eas. (l40> . Then, the dispersions for 
q{t),p{t) are proportional to 

(cosh(5±sinh(5sin(2tJo< + 6')), (41) 

distinguishable only by the + and — sign and constant prefactors, respectively. Therefore after the trans- 
port, the dispersions oscillate with twice the harmonic frequency and a relative phase shift of tt. The 
strength of this squeezing oscillation is thus solely ruled by the classical quantity 6. 

Finally, our classical initial conditions p{—to) — 1 /y/oJo together with the choice of the free parameters 
ci — C3 — LLjo/m,C2 — which we used to define the annihilation and creation operators in Eas. ( 1341331 . 
provide the correct initial dispersions of the quantum formulation in this approach 

{a\{Aq{^to)f\a) = xl/2 {a\{Ap{~to)f\a) = hy2xl. (42) 



6 Transport dynamics in a well-controlled regime 

In the following we consider an idealized situation for the transport, i.e. we assume that we could pro- 
duce arbitrarily shaped external potentials in the experiment while locally maintaining parabolic potentials 
around qo, i.e. 1 ^ \a[.^^{qQ)\/ujQ and |qo| ^ |ares(Qo)| for all positions qo or times qo{t). Deviations from 
these ideal conditions due to constraints in realistic trap configurations will be evaluated in Sect. In the 
ideal case we find from Ea.(l23> 

E{ti) = -e"^"'' / dt' e-"^°* iioit'). (43) 

J -to 

For to,ti 00 we arrive at the well-known result that the transferred energy corresponds to the squared 
modulus of the Fourier transform of the time-dependent force at frequency t^o 1 34 1 . Since we can decom- 
pose any function into a sum of symmetric and anti-symmetric parts go (^) = {(lo{t)+qo{^t))/2 + {qo{t)~ 
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go(-t))/2 = q^{t) +q^(t) we can write 



dt' sm{uJot')qo{t') 

to 



tl 

dt' cos{LUot')qo{t') 

'to 



(44) 



The two parts increase the amount of the transfeiTed energy independently. For a real transport, where start 
and stop positions differ from each other, we need anti-symmetric parts in the transport function. A simple 
conclusion from this is that any symmetric part of the transport function can only increase the transfeiTed 
energy while not contributing to the purpose of the transport, therefore we only need to consider anti- 
symmetric functions as candidates for transport, i.e. we take (0 = 0. By partially integrating Ea. (l43t 
two times and using initial conditions for the start and stop position and velocities Eq.(|2j, we can also 
rewrite the integral in Ea. (l43> as a direct functional of qo{t) that has a similar appearance but with an 
additional term. By symbols in this text with an extra tilde we denote quantities that are divided by the half 
of the transport distance &/2, e.g. qo{±to) = go(±^o)/(^/2) =±1. 

The average number of vibrational quanta transferred during transport can now be calculated from 
Eq.([l9j 

7(to) = mb^LOo\E{to)/LUo\y8h. (45) 

The energy increase in a transport therefore scales quadratically with the transport distance if the time 
span is fixed. Before we systematically study expression Ea. (l43t we will consider two examples for which 
analytical solutions exist. 

6.1 Two analytical examples 
First we take a sine function for the transport function qo{t) as used in the experiments described in II18I 

qo{t) ~ sin(t7r/2to) for — to < s < io- 

A graph of this function is given in FigQJ). Inserting q^it) into Ea.(l43> we find 

~ ~ 2 cosfxl 

S(to) = S(a;/tJo) = t^o — x phase with a; = cjo^o, (46) 

1 — (2a;/7r)" 

where we converted the time variables to dimensionless units, x = ojQtQ, so that x/2tt corresponds to the 
number of oscillation cycles. In these variables, |S(to)/tL'oP is independent of the frequency and plotted in 
Fig[D')- The energy transfer is decaying overall, but shows some oscillations arising from the dependence 
on the energy transfer on the phase of the internal oscillation at Iq. From Ea.(l45> we see that for an extreme 
nonadiabatic transport, i.e. a; ~ 0, we have gained the full potential energy of muj^lP' /2. Depending on the 
exact transport duration we observe regular intervals where the energy drops to zero and no energy remains 
in the internal oscillator's motion after the interval length 2<o- This is due to the phase sensitivity of the 
transport. From Ea.( l46> we have the proportionality |Sp cx (1 + cos(2a;)), so that we expect the first zero 
for x/2'iT — 1/4. However, for a transport in a harmonic well we need at least half an oscillation period for 
the ion to move to the other turning point, therefore we have instead x/2-k = 3/4 which is seen in Fig^) 
as the first root of the adiabatic factor The denominator in Ea.( l46t cancels the first root. The adiabatic 
energy transfer corresponds to the envelope of this function and is given by 4/(1 - (2x/7r)^)2. As we will 
see in the following the decay of the envelope can be sped up for different choices of the transport function. 

Typically, we want to have m6^a;o|S/a;oP/8ft < 1 in order to limit the maximum transferred energy 
to a few vibrational quanta. Let us consider some typical parameters for traps currently in use; we choose 
the axial frequency = 27r • 3 MHz, a typical average transport distance of about four traps (=control 
electrode widths), i.e. 6 = 400 //to, and to = 9.01218 u equal to the mass of a Beryllium ion. Then the 
adiabatic suppression factor should obey 

|S/c^oP<2-10-«. (47) 
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Fig. 1 Sinusoidal transport. In a) a normalized transport function for a sine transport is shown with a transport in 
2to — 3 time units. Figure b) shows the energy transfer as a function of 2;/27r = t/T, with T = 2'k/ijJ the period of 
the associated oscillation frequency. Zeros occur at the positions x/2n — (2n + 3) /4 with n = 0, 1, . . .. See text for 
details. 



In Fig[TJ)) we have not plotted the whole range until this criterion is fulfilled. It is satisfied for about 
x/2Tr > 30. Thus, for the given case adiabatic transport happens on a rather long time scale, i.e. durations 
of 2 • x/2tt = 2to/T = 60 cycles. The transport in the experiment | i8| which has used this transport 
function was performed over three times this distance requmng that |S/a;oP is lower by a factor of 9 
more. The adiabatic envelope has decayed to this value at about x/2tt ^ 52 yielding a transport duration 
of 2to/T = 104 cycles. Using a sine transport function the experimentally measured limit was around 
2tQ/T — 157 oscillation cycles (where 2to = 54 fis and ojq = Stt 2.9 MHz). This appears reasonable 
because the electrode array that was used in 1181 was rather sparse, thus not allowing for full control 
and maintaining the conditions assumed in this paragraph properly. Also, the envelope in this region is 
quite flat; so within the uncertainties of the experiment, the experimentally observed limit is in reasonable 
agreement with our estimation. 

We will look at an error function transport which turns out to be advantageous to the sine function in 
the second example 

qo{s) = Erf(2s/tp)/Erf(2to/ip) for - to < s < to, (48) 

where we renormalized it to arrive at the times ±to at the start and end position. In addition we have 
introduced another time tp which is nearly reciprocal to the slope of the transport function at the central 
point t = 0. Figl^t) is a graph of this function for {tp = 1, = 3/2) in arbitrary time units. Since we 
truncate the error function, we violate the second initial condition in Eq.(|3 in a strict sense. However, we 
are interested only in settings where tp/2to < 1, so that this constraint for the velocity can be satisfied 
arbitrarily well. The adiabatic suppression factor can be evaluated analytically 

E{x/ujq)/u;o - 2e-y^/^^Re{Eif[2x/y + iy/i]} /ETf[2x/y] x phase factor (49) 
^ 2e-y X phase factor, 

neglecting the part resulting from the finite initial and final velocities, and using the dimensionless variable 
y = LdQtp. Fig|2j)) illustrates the situation for y = 12 and x/2tt in a range of time intervals the same as for 
the sine transport but also satisfying 2to/tp > 1. It is clear that by using the error function the transport 
can be performed much faster than with a sine transport function, while still satisfying inequality Ea.( l47> . 
The full transport can now be performed in 2a;/27r ~ 6 cycles with tolerable energy transfer. Interestingly, 
taking the limit for large 2x/y ratios in Ea. (l49> removes the phase-sensitivity completely. However, we 
also want to note that the differences observed in these examples depend on experimental circumstances. 
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Fig. 2 Error function transport, a) graph of the renormalized error function, b) Adiabatic suppression factor for the 
error function transport. Details are given in the text. 



e.g. for very short transport distances, the adiabatic suppression factor does not have to be small. In 
this case the differences between the adiabatic suppression factors is marginal in a qualitative sense. This 
can be seen in comparing Fig[23) and Fig|2j5) for cases when only about |S/woP < 10"'^ is required, e.g. 
occuring for transport distances much less than an electrode width. On the other hand we find an interesting 
and advantageous distance scaling behaviour from Ea.( l49l i: transporting longer distances does not require 
much longer time intervals. For example, Steane |20| estimated that within large-scale operation for the 
processing of a typical gate an average transport distance of ^ 22 traps is needed. By employing an error 
function transport we find that this is feasible with less than a quantum of transferred energy using the 
parameters (y = 13,a;/27r = 4), i.e. in already about 8 oscillation cycles, only about a third more time 
than for a transport over 4 traps. The average velocity for such a transport is then considerably higher. 



6.2 Near-optimum transport functions 

In an attempt to optimize the transport function we can expand expression Eq.(l49> up to the first order 
correction L31J 



E{x/uJo)y ~°° 2u>o ( e ——=== cos(2r& + arctan(6/r)) ) /Erf [r] x phase factor, 

(50) 



with r = 2x/y = 2to/tp, b = y/A. Because of 2r6 — a; the zeros of the suppression factor |Sp are equally 
spaced as in the sinusoidal example if we disregard the phase in the cosine function, i.e. half periods 
A(x/27r) = 1/2 since the particle can arrive from two different turning points at the end of the transport. 
From this expansion it is clear that the ratio r basically determines the magnitude of the second term on the 
rhs of Ea.( l50> and suppresses the phase sensitivity as it increases. If r is chosen large enough the energy 
transfer is dominated by b. To find some conditions that are close to optimum we proceed in the following 
way: first, from the transport distance and achievable frequency we can evaluate the upper bound for the 
adiabatic suppression factor as in Ea. (l47> . Because we have to satisfy 2to > tp we then choose b large 
enough to suppress the first exponential factor to fulfill the given criterion. This procedure defines the 
asymptotic value of energy transfer for large x/27r as shown in Figl^Ji). We then choose the interval length 
x/27r = rb/iT by defining r such that we are just in the asymptotic range. The near phase-insensitivity can 
then be thought of as a result of the extremely slow start where the phase information in the limiting case 
r ^ cx) in Eq.(l49> gets totally lost. 
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6.3 High-frequency limit, adiabatic transport, and approximate trajectories 

To better understand the behaviour discussed in previous sections, we present a few more general consid- 
erations. We write the time-dependence of the transport function according to qo{t) — ^o{t/to) = ^o{t) 
so that i9o only depends on the dimensionless variable r (for the error function example we also keep the 
ratio tp/to fixed). Making the substitution t' ^ r in the integral in Ea.( l43> we find 

E{to)/ujo = {oJohy'^ J dr e-''^°^'>-^do{T) x phase factor (51) 

wotoj>oo 2 ^ ^ ^:—4t cosfwoio + n7r/2) x phase factor. 

The exponent in the integral relates the two time scales in ujq and to. Using the method outlined in appendix 
|A]we expand it into the sum given in the second line in the limit ojo^o oo assuming that is sufficiently 
smooth. In this expansion the derivatives at the start position (and end position due to anti-symmetry) define 
the energy transfer in the transport, and thus fully characterize the transport function for the transferred 
energy in the adiabatic limit. This provides us with a reason for the difference we observed above for the 
error function and sine examples. The second derivative for the sinusoidal transport is nonzero at ztto and 
much larger than in the case of the error function. In the latter all derivatives are damped by a gaussian while 
the ones for the sine transport alternate. Furthermore, we see that we can in general decrease the transferred 
energy for larger values of the product ujoto, i.e. by taking ojq to inifinity (high-frequency limit), we can 
lower the adiabatic suppression arbitrarily, on the other hand, slowing down the motion by increasing the 
length of the duration of the transport 2to, we move into the adiabatic regime. For infinitely slow motion 
we end up with zero transferred energy. These two limiting cases are formally equivalent because the 
energy transfer depends only on their relative time scale. We can perform the same expansion starting 
from Ea.( l43t directly and use the relation Uc{t) = — E* {t)]/2iuJo to find approximate trajectories 

valid in the same limits 

= - E -1+2 
n=0 ^0 



COS (^) q^-^'\t) - COS + to) - ^) q^'H-to) 



7 Regularized trap-electrode waveforms, potential fluctuations and aspect- 
ratio rule 

7. 1 Determination of waveforms 

So far, we have said nothing about how to determine the waveforms applied to the electrodes. As soon as 
we have the waveforms at hand for a given model electrode configuration we can determine the magnitudes 
of perturbations. This is done in the next section. Here, we seek optimum solutions for a given electrode 
structure in order to keep the uncontrolled part (f>res{q) of the total potential in Eq.@ small. The time- 
dependent electric potential is created by a linear superposition of the available control potentials <j)m (<?) 
and dimensionless time-dependent amplitudes (t) of the form 

t) = E a„i{t)(p,n{q)- (52) 

m 

To optimize waveforms for the time-dependent amplitudes for the transport problem we find a measure of 
the discrepancy by integrating over the residual non-matched part according to 

rqo{t)+Sq 

min / \4'ies{q,t)\'^ dq V t (53) 
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while (j){q, t) from Ea. (l52> enters here through Eq.(|4}- For any time t we want to find a set a„i, cpo for which 
expression Ea.( l53t is minimal. The integration is performed over an interval moving with the minimum 
of the parabolic potential well, i.e. [qQ{t) — Sq,qo{t) + Sq] and assuming a unity weight factor in the 
integrand. We do not consider in this range any lag of the ion due to acceleration and deceleration since for 
an adiabatic transport and experimental conditions the lag is much smaller compared to the optimization 
range, ipo (t) represents here another degree of freedom that does not perturb the dynamics but might allow 
one to more optimally choose the harmonic potential well by arbitrarily offsetting the desired parabolic 
potential for best fit. Condition Ea.( l53> is readily converted into a linear system of equations by taking 
partial derivatives for the amplitudes a„i and (po, and setting them all equal to zero. The minimization 
problem in Eg.dSSt then reads 

^ amit) J dq (t>^(f>j +ipo{t) J dq ^ j dq [q ~ qo{t)f(l)j V (54) 

Sa S K 

where we dropped the explicit integral bounds and arguments of the potentials (pmiq) for the sake of 
simplicity. Bold symbols denote matrices, underlined symbols vectors. 

The optimization problem can then be formulated in terms of the linear system 

S+o • a+o ^ ( Sa ) • ) - VK \ft, (55) 

with r] = mti>^/2Q. All quantities are functionals of qo and for a given transport function qo{t), we need 
to solve the equations at every point in time. As a result we obtain the waveforms a„i (t) = a„i [qa (t)]. 

Typically we choose an optimization range of 2Sq ~ 0.5W^ for electrodes of width W. This is usually 
much smaller than the mean distance between most of the contributing electrodes to the center of the 
parabolic well. Thus, due to the slow decay of the axial potentials, the curvatures of distant electrodes are 
similar, and their contribution differs locally only by a multiplication constant. This is particularly true 
for experimental situations where the high electrode density typically makes the system Ea.(l55> nearly 
singular. A straightforward least-square method, such as 

a+o = argmin{||S+o -a+o - TyXll^} , (56) 

is therefore not well suited for finding waveform amplitudes. For high electrode density, a tiny step 
qo + S might change the individual electrode amplitudes exponentially fast. In these cases the matrix S+g in 
Eq.(|55} filters out too much information from a+Q to invert this system properly. In mathematical language 
these kind of problems belong to the family of discrete ill-posed problems that can be numerically solved 
using regularization approaches |35|. Here, the lost information is fed back in the minimization process 
via a Lagrangian multiplier concerning the smoothness of a„i{t), or curvature etc. in amplitude space. If 
we apply a Tikhonov regularization to the given problem we have to solve 

argmin{||S+o - a+o -'/^IP + '^^l|L(a+o -a;o)|P} (57) 



in order to determine smooth time-dependent waveform amplitudes am{t)- In Ea. iSH the regularization 
parameter i' corresponds to a weight factor between the original least-square minimization and the addi- 
tional side constraints, while a+g can be used to find solutions near a prescribed setting. The smoothing 
properties of this optimization originate from a common and simultaneous minimization of both terms. 
L is a linear operator that can be used to feed back different kinds of information to the amplitudes. For 
the results given here we took for L the unity operator, and also a+p = 0. Since we only want to limit 
the amplitudes a„i to some appropriate experimental values and stabilize the solution, our interest is not 
to determine the overall minimum of Eq.(|57} in a self-consistent way. For our convenience we choose v 
manually to make the parameters compatible with available technology. 

We can summarize the advantages of these methods to the current optimization problem: 
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1 . The regularization method selects only nearby electrodes for creating a local parabolic potential, and 
disregards tiny linear contributions from distant electrodes which would require large amplitudes to 
effect small changes. 

2. The choice of the regulariyation parameter Umits the amplitudes a,n to practical experimental values. 

3. It is robust against changing the electrode density (here, governed by the widths W). This will be of 
importance in the next section. 

4. It stabilizes the output waveforms and smoothes sharp features in the time-dependence of the am- 
plitudes. Different constraints can be set via the L operator, defining bounds or curvatures in the 
amplitude space. 

For more detailed information we refer the reader to the mathematical literature 1351 . A typical example 
of a parabolic potential created through superposition of an array of electrodes and the time-dependence of 
amplitudes is shown in Figl^and further discussed in the next section. 



7.2 Potential fluctuations and aspect-ratio rule 

Based on a reasonable multi-electrode structure we want to estimate how well we can meet the require- 
ments on transport potentials stated above, in particular, how stringently we can meet \a[.^^{qo)\ / luq <C 1 
and I Ores (go) I *C l^ol- We employ the definition of waveforms and the method from the previous section 
for extracting the residual, uncontrolled potential 0res from which we perturbatively derive the effect of 
imperfections on the transport. As a simple model electrode structure for transport in single and multi- 
layer traps, we use the "railway track" electrode configuration sketched in Fig|3t) which might be a simple 
model for transport in single and multi-layer traps 1191 1221 l32l 1331 . The transport occurs along the long 
arrow where we assume the ion is held radially by RF fields and controlled axially by the electrical fields 
arising from the potentials of the "stripe" electrodes depicted in Fig|3^). We are mainly interested in the 
scaling behaviour as a guideline for general design rules. Waveforms that are actually used in experiments 
should be based on more accurate numerical potentials and generalized versions of Ea.(l54> for all three 
dimensions. We are finally interested in the trade off between adding electrodes, by shrinking the electrode 
distances/widths along x, and the amount of control that is gained in that way. 

We can model this arrangement as a sum over the potentials <j>rn (x) of several infinitely long (in the 
y-direction) stripe electrodes that are distiibuted along the x-axis 

ElV ( W \ 
am0m(i) = — > a,„ arctan -j— -j— — ■ , (58) 

where each (/)„i [x) is the exact solution of the Poisson equation for an infinitely long stripe at position m • W 
that is embedded in a ground plane. For convenience we choose for the individual potentials in this basis 
set a potential on the electrodes of IVolt — \ V. We denote symbols with a hat as quantities normalized 
to the ion distance Zion to the surface, e.g. the normalized electrode width W — W/ Zioa- Figjsj)) shows 
the behaviour of {x) for various geometric aspect ratios W . We see that a plateau-like structure starts to 
form for ~ 2 and larger resulting in small field gradients along the transport direction in the center of 
each stripe electrode. The maximum frequency at the center of electrode m = is obtained from 



,2 



Q d' 

Zion) = 2~^T^<^o(i) 



ion 



-g^li^,,,.,^,. (59) 



The proportionality uj cx is rather unusual and stems from the fact that we scale only a single dimension 
(along x). The last factor of the second equation is solely defined by the aspect ratio W and thus by the 
geometry of the trap. It exhibits a maximum for ^ 0.78 and decreases only significantly for small width- 
distance ratios W < 0.5. Using the mass of Beryllium, W = 1, amplitude oq = —2 and Zio,, = 40 fim as 



Copyright line will be provided by the publifiher 



16 



R. Reichle et al.: Transport dynamics of single ions in segmented microstructured Paul trap arrays 




Fig. 3 a) The control electrode array model discussed in the text with x the axial coordinate, Zion the vertical distance 
of trapped ions to the electrode plane, and W width of an electrode. No RF electrodes are shown to hold the ion in 
radial directions. The electrode 'stripes' are assumed infinitely extended along ±j/, and arbitrarily many along ±x. 
Ions are transported along the long arrow from —6/2 to +6/2. b) Modulus of the electric potential of a single electrode 
at height Zion along x (relative to the center of that electrode and normalized to Z[on). The electrode is embedded in an 
infinite ground plane and held at — 1 1^, and W = W/z[on- 

the ion surface distance for the surface-electrode trap as used in f32], we find an axial frequency of about 
Wo ~ Stt • 9 MHz. We chose this low value for |ao| to be compatible with typical maximum voltages as 
created by CMOS electronics 1191 . With these parameters we have created waveforms for the transport 
of a confining harmonic well utilizing the regularization approach of the last section. We used a set of 41 
electrodes while the transport was over four electrode widths, b ~ A, around the central electrode 21 of 
this array. For the transport we used the error function Ea.(l48> of Sect|6lwith a transport duration of 8 
oscillation cycles x/2Tr = 4 and y = 12. We then determined the lowest order deviations from an ideal 
harmonic potential with constant trap frequency and controlled acceleration, uj{t)/ujo and —aye^{qo)/ujQ, 
respectively, in Eq.(|6|l for the aspect-ratios W — (0.5, 1.0, 1.5, 2.0). Figs|^) and b) illustrate an example 
for the superposed potentials, and for a set of waveform amplitudes for an error function transport in 8 
oscillation cycles. 

The choice of the regularization parameter i/ is not obvious, because we have to deal with a set of near- 
singular matrices all at once. As mentioned earlier we do not aim for self-consistent methods to determine 
1/ and an absolute minimum of the expression Ea.( l57> [35J. In our context we are more interested in a 
feasible implementation compatible with given experimental constraints. The choice of the regularization 
parameter affects both the stability of the linear system and the size and smoothness of the amplitude vector 
Qm ■ In a strongly regularized inversion more stability is added to the solution, forcing the amplitudes to 
be of limited size. Because of this bound the solution can not closely approximate the desired shape of the 
potential anymore, so the deviations from the ideal case increase. A weak regularization scheme, on the 
other hand, adapts more closely to the desired potential shape, but reveals random fluctuations and noise 
on the solution waveforms a,n{t). Also the singular behaviour increases dramatically with an increase 
of the number of electrodes, and larger parameters ly have to be chosen. This latter property makes a 
direct comparison of the results among various aspect-ratios W difficult. Nevertheless, we can make some 
qualitative and general statements. 

The results for our sample configuration are summarized in FiglS] The upper graph in FiglSjdisplays the 
uncontrolled acceleration —ai-es{qo)/ujQ, and the middle and lower graphs the frequency modulation uj{t) 
for various aspect ratios. In both figures a dramatic change of the curves is observed around W ^ 1. While 
for smaller ratios the frequency fluctuations are in the percentage range (middle panel), the emulation of 
the potential for larger ratios is much worse due to the constraint |a„i| < 2 (lower panel). Frequencies 
drop by more than 90% already for the W = 2 calculation. Only in the strongly regularized scheme, 
did we find a direct correlation of the solution vector a„i to the choice of the regularization parameter. 
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x/W 

~-4 -2 2 4 

t/T 

Fig. 4 a) Creation of a harmonic well by superposing potentials of the electrodes of the stripe configuration of Fig|3ji)- 
The fine lines represent individual electrode potentials that sum up to the total potential represented by the thicker line. 
The dashed-dotted line shows the ideal harmonic potential and the two vertical lines indicate the range of optimization 
at this given location. By symmetry the extrema of the individual electrode potentials are located at integer values of 
xjW , with the maximum of electrode potential 21 at the origin, b) Waveform amplitudes of electrodes 18 to 24 for a 
transport from trap 19 to 23 for W — 1. Waveforms from more distant electrodes still contribute but are not shown for 
the sake of clarity. Amplitudes 21-24 are marked by symbols for better visualizing their traces. The transport is done 
according to an error function, while the abscissa represents time in units of oscillation periods t/T. 



In the weakly regularized scheme the amplitudes were limited by other lower bounds, and the waveform 
solution appeared similar over a large range of i^, but exhibited a much noisier behaviour. This enhanced 
sensitivity is an indication that inclusion of more electrodes (smaller W) does not improve the quality of 
the solution anymore. The linear system becomes more singular and exhibits more rank-deficiency, i.e. 
rows and columns become more equal and their inclusion adds more redundancy. For the given parameters 
we observed that for W — 0.5 the transition from a regularized to a weakly regularized solution occured. 
Therefore, our results indicate that W = 0.5 — 1.0 should be optimal for the configuration discussed here. 
For larger aspect- ratios we found that the coverage of curvatures of the individual potentials along the 
transport axis is not sufficient for the necessary amount of control. 

The other constraint, i.e. |ares(9o)| \qo\, of a controlled transport force, has to be interpreted dynam- 
ically. Since the acceleration force depends on the time duration in which the transport is performed, this 
requirement can be violated for a slower transport. In Fig|3 we show that in the initial phase the pertur- 
bations overwhelm the transporting acceleration for aspect ratios W > 1.5 or larger. Results for smaller 
aspect ratios are not given in this figure because the transport force by far dominates the excess force and 
lead to a fully controlled transport. 

In general, fluctuations in the frequency and transport force affect an energy transfer according to 
Ea.( l23> . This introduces violations to the symmetry of the transporting force and leads to an enhanced 
energy transfer as seen from Eq. J44> . We have not included higher order terms in our discussion, because 
we aim for experimental conditions to perform a transport in the well-controlled regime. However, they 
are inevitable for longer transport distances and other types of motion, such as nonadiabatic transport, or 
splitting of ion groups where they might lead to large energy transfers. 



Copyright line will be provided by the publifiher 



18 



R. Reichle et al.: Transport dynamics of single ions in segmented microstructured Paul trap arrays 



CVJ o 


-4 


3 


0.05- 




0- 






0) 

<ns" 
1 


-0.05- 




1.02r 




1 01 . 




o 0-8 
3 

^ 0.6 
3 0.4 
0.2 




t/T 



Fig. 5 First and second order perturbation to the transport (thiclc lines), and their dependence on various aspect 
ratios W. The abscissa t/T corresponds to the number of oscillation cycles with period T = 2tt/uio. The labels on 
the graphs denote the value of W. (Upper panel) Residual, uncontrolled acceleration — ares(go)/i^o (dimensionless) 
which should be compared to (jo/i^o (dotted line). This is only given for the large aspect ratios W = 1.5, 2.0. The 
graphs for smaller ratios would be close to zero on this scale. (Middle and lower panels) Relative frequency modulation 
uj{qo)/uJo during an error function transport for various aspect ratios. For explanations see text. 



8 Conclusions 

In conclusion, we have analyzed the dynamics of single ion transport in microstructured linear Paul trap 
arrays. We have modeled the transport by a forced and parametrically excited harmonic oscillator and 
have presented a theoretical framework for its description. We have derived exact analytical expressions 
for the classical as well as quantum dynamics and reviewed their related properties. In particular we have 
expressed the Heisenberg operators by the approach of Kim et al. [21 J through the dynamical quantities of 
the related classical solution. We have given explicit analytical expressions for the classical energy transfer 
involved in these transport phenomena and derived expressions for the lowest order deviations from ideal 
transport that will necessarily appear for unfavourable ion trap layouts. For current trap technology we 
have evaluated durations for a fast adiabatic transport and found that they depend strongly on the external 
force employed in the transport. According to these results, the adiabatic single ion transports of reference 
[18,1 could be sped up by more than an order of magnitude with negligible energy transer to the motion. We 
determined appropriate transport waveforms and found that with an adiabatic transport over four electrode 
stripes of size roughly equal to the distance of the ion to the nearest electrode and frequencies in the range 
of ~ 9 MHz is feasible in about 6 oscillation cycles. Our results also indicate that a full control over the 
transport is available, where perturbations to a harmonic oscillator potential are negligible at aU positions 
and times. By directly relating deviations from these ideal potentials to the aspect ratio of the trap, we 
have found a practical design rule that should be valid for trap layouts more general than the one given 
here. The ratio of a control electrode width to the distance to the ion should be in the range 0.5 — 1 
for a well-controlled regime. Our example suggests that a higher electrode density does not appreciably 
improve transport performances. This provides important insight into the amount of resources needed to 
realize large scale implementations of ion trap based quantum computers. Transport in a confining well 
of constant frequency might also enable continuous cooling processes during the transport. If eventually 
experiments allow one to maintain a well-controlled regime during the transport, performing quantum 
processing during transport is conceivable, possibly leading to appreciably shorter processing times. 
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A Integral expansion 

We employ a mathematical theorem proven within the formalism of h-transforms, see for example theorem 
3.2 of |36|: If .g(T) has 1 continuous derivatives while (7(^+2) is piecewise continuous on the real axis 
[a, b] then 



/(A) 



n=0 ^ 



iXa 



(60) 



If we also require lim^^a+ = limr^fc- .9^"H'^) = for n = 0, .., iV - 1 it holds that /(A) 

o(A-^). 
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